R2 = 1000;
R0 = 1000;



simulacion_10 = dlmread('-10V.txt','\t',1,0);
simulacion0 = dlmread('0V.txt','\t',1,0);
simulacion2 = dlmread('2V.txt','\t',1,0);

figure
hold on

simulacion2(:,2) = simulacion2(:,2) - simulacion2(1,2);
simulacion0(:,2) = simulacion0(:,2) - simulacion0(1,2);
simulacion_10(:,2) = simulacion_10(:,2) - simulacion_10(1,2);

tau2 = calcular_tau(simulacion2, 2, 4)
tau0 = calcular_tau(simulacion0, 2, 4)
tau_10 = calcular_tau(simulacion_10, 2, 3)

plot(simulacion2(:,1),simulacion2(:,2),'g-','Linewidth',2)
plot(simulacion0(:,1),simulacion0(:,2),'r-','Linewidth',2)
plot(simulacion_10(:,1),simulacion_10(:,2),'b-','Linewidth',2)

legend(sprintf('simulacion con Vmin = 2 V, Tau = %e seg', tau2),
		sprintf('simulacion con Vmin = 0 V, Tau = %e seg', tau0), 
		sprintf('simulacion con Vmin = -10 V, Tau = %e seg', tau_10),
	    'Location','northoutside')

xlabel('Tiempo [seg.]')
ylabel('Variación de tensión en el diodo [V]')

grid minor

print('grafico.png','-dpng');
